Many-body theory of surface-enhanced Raman scattering 
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A many-body Green's function approach to the microscopic theory of surface-enhanced Raman 
scattering is presented. Interaction effects between a general molecular system and a spatially 
anisotropic metal particle supporting plasmon excitations in the presence of an external radiation 
field are systematically included through many-body perturbation theory. Reduction of the exact 
effects of molecular-electronic correlation to the level of Hartree-Fock mean-field theory is made for 
practical initial implementation, while description of collective oscillations of conduction electrons 
in the metal is reduced to that of a classical plasma density; extension of the former to a Kohn- 
Sham density-functional or second-order M0ller-Plesset perturbation theory is discussed; further 
specialization of the latter to the random-phase approximation allows for several salient features 
of the formalism to be highlighted without need for numerical computation. Scattering and linear- 
response properties of the coupled system subjected to an external perturbing electric field in the 
electric-dipole interaction approximation are investigated. Both damping and finite-lifetime effects 
of molecular-electronic excitations as well as the characteristic fourth-power enhancement of the 
molecular Raman scattering intensity are elucidated from first principles. It is demonstrated that 
the presented theory reduces to previous models of surface-enhanced Raman scattering and leads 
naturally to a semiclassical picture of the response of a quantum-mechanical molecular system inter- 
acting with a spatially anisotropic classical metal particle with electronic polarization approximated 
by a discretized collection of electric dipoles. 

PACS numbers: 71.10.-w, 33.20.Fb, 31.15.xp 



I. INTRODUCTION 

Inherently exceedingly weak, Raman scattering [Ij] of 
incident electromagnetic radiation from a molecular tar- 
get occurs, approximately, only once out of every million 
photon-molecule scattering events. These few inelasti- 
cally scattered photons carry away a fraction of energy 
less or more than they had originally with the difference 
being deposited into or liberated from molecular vibra- 
tional, or, to a lesser extent rotational, excitation. Resul- 
tant changes in electronic polarizability with respect to 
underlying nuclear geometry ultimately are encoded in 
spectral fingerprints of molecular structure. Discerning 
such small Raman signals from a large elastic Rayleigh- 
scattering background is challenging, yet potentially re- 
warding, as Raman spectra can reveal detailed molecu- 
lar structural information unresolved by other spectro- 
scopies. 

In the 1970s it was first observed [2| and later un- 
derstood [3, [1] that certain molecules, if adsorbed onto 
roughened noble-metal substrates having structure on 
the subwavelength scale, experience a surface-enhanced 
Raman scattering (SERS) of incident photons that man- 
ifests itself in approximately a millionfold boost in sig- 
nal in comparison to the normal Raman effect. Out of 
this new phenomenon, the field of surface-enhanced Ra- 
man spectroscopy was rapidly born; early reviews can be 
found in Refs. @, @, 0j @]- Two mechanisms are gener- 
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ally believed to account for enhancement of the Raman- 
scattered field: one of chemical and the other of electro- 
magnetic origin. The former rests upon the idea that a 
chemical bond is formed between the adsorbed molecule 
and the metal, thus allowing for charge-transfer exci- 
tations to occur between the two systems. The latter, 
which is widely regarded as being dominant involves 
the coupling of external radiation to surface-plasmon ex- 
citations at the metal-dielectric interface; plasmons are 
quantized collective oscillations of conduction electrons 
against the positive ionic background that can act to 
enhance and focus incident light to subwavelength di- 
mension below the diffraction limit. When optically ex- 
cited, these plasmons broadcast their enhanced field to 
nearby Raman-active molecules which inelastically Ra- 
man scatter photons back to the metal. The scattered 
photons can recouple into the plasmon modes of the 
metal and, subsequently, be rebroadcasted toward a de- 
tector. Both normal and surface-enhanced Raman scat- 
tering events are linear processes (depending only upon 
the linear polarizability) yet, if both incident and Raman 
scattered frequencies are resonant with plasmon excita- 
tions in the metal, their enhancements can multiply to- 
gether to yield a fourth-power enhancement of the normal 
Raman-scattering intensity. 

Recently, surface-enhanced Raman spectroscopy has 
gone through a renaissance which is, in part, attributable 
to the first observation of single-molecule SERS in the 
late 1990s [Tol . [T]| . Exhibiting a giant boost in Raman 
signal (by 10 or more orders of magnitude in comparison 
to the normal Raman scattering from a single molecule 
in free space), single- molecule surface-enhanced Raman 
spectroscopy provides an even richer variety of molecular 
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structural information that is free from ensemble averag- 
ing. With the ability to measure single-molecule Raman 
spectra, basic research has continued to progress under 
the impetus of utilizing SERS techniques, inter alia, as 
an ultra sensitive analytical probe having broad utility 
in the biological sciences; see, e.g., Refs. ITsl. IT3| . 

Today, over thirty years after its initial discovery and 
a decade after the observation of single-molecule SERS, 
significant work is underway to systematically character- 
ize the best conditions for single-molecule SERS activity 
in individual and arrays of nanoscale metal particles as 
a function of size, shape, and interparticle separation, 
among others, with respect to the wavelength of light 
[TH - Independent variation of each of these nanoscale 
characteristics is controllable in the laboratory, with the 
number of permutations exceedingly large. Guidance 
from predictive theory would be of immense utility in 
this pursuit, however, current theoretical methods can- 
not completely meet this challenge. 

Involving the coupling and interaction of a Raman- 
active molecular system with one or more nanoscale 
metal particles under the influence of an external radia- 
tion field, the SERS effect presents a complicated many- 
body problem blending together concepts from quantum 
chemistry and molecular spectroscopy, condensed-matter 
physics, and electromagnetism. Undoubtedly, its present 
incomplete theoretical description is rooted in the com- 
plexity of its basic processes. As it not our intent to ex- 
haustively summarize thirtyyears of theoretical progress, 
we defer to the reviews d, H S H, 0] and briefly discuss 
only a few recent and notable approaches based upon the 
complementary starting points of classical and quantum 
mechanics [16|. Almost all previous approaches can be 
placed into one of these two categories. 

The optical and plasmonic properties of spatially 
anisotropic nanoscale metal particles (and particle ar- 
rays), which may have dimensions from tens to hundreds 
of nanometers, are well described within classical elec- 
tromagnetic theory. Useful physical information can be 
gleaned from electrostatic model calculations [13, [III as 
well as from full numerical solution of Maxwell's equa- 
tions [l^, [131 , such as the magnitude and location of en- 
hanced electromagnetic fields located near the particle's 
surface [2lj: so called electromagnetic hot spots. How- 
ever, while a classical description based upon the metal's 
underlying continuum dielectric function may be appro- 
priate for a nanoscale particle, an adsorbed molecular 
system undergoing inelastic Raman scattering of photons 
is properly described only from a microscopic point of 
view. 

To this end, recently, a fully quantum description 
of the combined nanoparticlc-molecule system has been 
developed within the Kohn-Sham framework of time- 
dependent density- functional theory [13]. Both particle 
and molecule are represented by the same total wave 
function. Leaving aside shortcomings inherent in the 
choice of electronic-correlation functional and its non- 
systematic improvability, this approach is limited, due 



to computational restrictions, to the treatment of small 
metal particles containing, at most, on the order of one 
hundred atoms [1^. With such small numbers, metallic 
particles display the discrete electronic-excitation struc- 
ture more typical of clusters than of the bulk-like res- 
onance continuum exhibited on the nanoscale. In this 
single-particle picture, molecular electronically-excited 
states have an infinite lifetime as there is no mechanism 
for their damping due to the presence of the metallic 
system; this is a consequence of the fact that there is no 
explicit treatment of the interaction between molecule 
and particle beyond that specified in the single-particle 
Kohn-Sham formalism. Ad hoc empirical parametriza- 
tion is employed to mimic these basic interactions, and, 
in turn, damp molecular-electronic excitation. In this 
way, deficiencies in the theory are corrected, leading 
to p redictions which compare sensibly with experiment 
[22l[23. [25I . Two additional notable quantum approaches 
to SERS based upon density-matrix calculations have 
recently appeared in the literature [2^, [13] ■ The inter- 
action of a Raman-active molecule supporting two elec- 
tronic states (plus several vibrational substates) with two 
nearby nanoscale Ag spheres, described through an ex- 
tended Mie theory, is presented in Ref. [2y]. With 
appropriate choice of parameters and inclusion of phe- 
nomenological damping mechanisms, Raman-scattering 
cross sections enhanced by 10 orders of magnitude are 
demonstrated in comparison to the normal Raman ef- 
fect. Second, in Ref. [13], enhanced resonance Raman 
phenomena are studied within a combined eight-state 
density-matrix approach where the molecular subsystem 
is represented in a four-state basis involving molecular 
ground and electronic, vibrational, and electronic and 
vibrational excited states while the particle subsystem is 
represented in a two-state basis consisting of ground and 
excited states. Within this model, it is predicted that 
the largest resonance Raman enhancements occur when 
a molecule, which absorbs light far from the particle's 
resonance maximum, is excited at the resonance maxi- 
mum of the particle. This prediction is believed to be 
caused by the shifting of molecular resonances due to the 
strong coupling between molecule and particle. 

Each of the approaches described above have both pos- 
itive and negative attributes. It is therefore natural to 
envision blending their best features and, simultaneously, 
to explicitly treat the coupling between molecule and par- 
ticle so as to avoid parametrization. It is the purpose of 
this article to do exactly that. 

Here we present a formal ab initio many-body theory 
underlying a unified and didactic approach to the de- 
scription of single-molecule SERS from a nanoscale metal 
particle at zero temperature. Emphasis is placed on de- 
veloping a rigorous, yet computationally tractable for- 
malism. In anticipation of practical initial numerical im- 
plementation, specialization is made, within the Born- 
Oppenheimer approximation, to a Hartree-Fock (HE) 
mean-field description of the electronic states of the 
molecule, while quantized collective oscillations of metal- 
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lie conduction electrons are described by their classical 
plasma density. These approximations, as we have ap- 
plied them within our minimal model, limit the possi- 
bility for charge transfer between molecule and metal, 
and effectively restrict our current presentation to the 
electromagnetic mechanism of SERS. We point out that 
a program similar to the minimal implementation of our 
approach has already been introduced in the time domain 
at the level of time-dependent HF theory for the molec- 
ular system and an electromagnetic boundary-element 
method for the particle [28|. Our approach is comple- 
mentary and more general in the sense that we develop 
the full many-body theory starting from the exact many- 
body Hamiltonian for the interacting molecule-particle 
system and its interaction with an external electric field 
in the dipole approximation. Using a generalized M0ller- 
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Plesset perturbation theory [29|, the effects of interac- 
tion with both particle and field are systematically and 
explicitly built into both nonperturbative and pertur- 
bative expressions for the molecular-electronic Green's 
function. Electronic-correlation effects beyond HF the- 
ory such as those of nth-order M0ller-Plesset perturba- 
tion theory may be rigorously and straightforwardly in- 
cluded, while, alternatively, it is also clear how to treat 
the molecular-electronic sector of our theory within a 
Kohn-Sham formalism (soj . Due to its general formu- 
lation, our approach recovers other approximate results 
from the literature, and, further, admits certain well- 
known observable features analytically upon invoking an 
analytic model for the particle's response. In particu- 
lar, a closed-form expression for the quantum many-body 
SERS intensity, which is of the generic form 



k'X(-?i^k')-(^Sl(Qj-Qo)-VQ,5;;^,,,,(fiL^k,-fi^k')kj)-9rg(?i^k)-4^^(k) 
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with polarizability transition moments 5j\/ , normal mode 
coordinates Qj, and incident and Raman-scattered en- 
hancement factors g and g' , is developed and presented 
below in Eq. ((79|) . To our knowledge, this is the first 
place in the literature where a SERS intensity is de- 
rived, entirely from first principles, that explicitly treats 
the coupling and back reaction of a quantum molecular- 
electronic system with a nearby metallic particle in the 
presence of external perturbing radiation. 

In Sec. [ni we review an early and insightful classical 
model of SERS based upon the couphng and interaction 
of two dipoles with each other and with the external elec- 
tric field. Two basic and essential types of interaction, 
the image effect and the local-field effect, are discussed 
and used to motivate our quantum-mechanical general- 
ization. Starting from the exact molecule-particle Hamil- 
tonian, the quantum many-body theory of SERS is devel- 
oped in Sec. Illll where the effects of interaction of molecu- 
lar electrons with metallic conduction electrons and with 
the external electric field are built into the underlying 
molecular-electronic Green's functions. The former ef- 
fect, which accounts for the repeated interaction of the 
molecule with its own image, is included to infinite order 
in Sec. IIII C H while the latter effect, which describes the 
interaction of the molecule with the local electric field of 
the particle, is included perturbatively in Sec. IIII C 21 
The random-phase approximation of the particle's polar- 
izability is invoked in Sec. IIII C 3l in order to demonstrate 
certain key properties analytically. Connection between 
the Green's function and scattering T-matrix is reviewed 
in Sec. IIII Dl and, subsequently, allows for the quantum- 
mechanical normal Raman-scattering intensity and en- 
hanced Raman-scattering intensity to be computed in 



Sees. IIII D II and IIII D 21 respectively. Equation 
which displays a first-principles quantum-mechanical ex- 
pression for the SERS intensity, is a major result of our 
work. Lastly, in Sec. IIIIE H linear-response theory is 
reviewed and used to compute the induced density of the 
interacting molecular system in Sec. IIIIE 21 and its in- 
fluence upon the dynamics of the conduction electrons of 
the particle in Sec. IIII E 31 Two appendices are devoted 
to the inclusion of molecular electron-electron interaction 
effects with density-functional theory, and to the Green's 
function based definition of the polarization propagator 
and linear polarizability. 

Summation is implied over all repeated Greek indices. 
All integrals of the form J (Px are taken over the vol- 
ume of all space, while those of the form J d'^x = J d^xdt 
are taken over the volume of all space and over all times 
from negative to positive infinity. Further, unless other- 
wise indicated, all time and frequency integrals run from 
negative to positive infinity. Molecular electronic state 
labels k,l, . . . refer to occupied or hole states, labels 
a, b,c,d, . . . refer to unoccupied or particle states, while 
the labels p, q,r^s, . . . are reserved for unspecified states. 



II. REVIEW OF CLASSICAL MODEL OF SERS 

Enhanced Raman scattering from a molecule can al- 
ready be described to some extent at a classical level of 
theory. In 1980, Gersten and Nitzan ^ proposed a sim- 
ple model consisting of two interacting electric dipoles: 
one dipole d'-'^' representing a molecule and the other 
p(i) representing an arbitrary polarizable body (taken 
here to be a metal particle) located nearby; see Fig. [T] 
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FIG. 1: (Color online) In the SERS process, an incident exter- 
nal electric field Eo induces first-order molecular and particle 
dipole moments d''^-' and p'^', which rebroadcast the incident 
radiation and, in doing so, couple to each other. Such a pic- 
ture, which includes the reaction and back reaction of the 
dipoles to the incident field and among themselves, provides 
a insightful classical description of SERS. 



Both dipoles are induced, at first order, by an external 
electric field Eq that varies harmonically in time with 
frequency and have corresponding molecular and par- 
ticle polarizabilities gm = <5a/[Q] and ap, where Q la- 
bels a nuclear vibrational coordinate of the molecule; Hm 
and ap may additionally depend upon frequency; tildes 
denote Fourier inversion to the frequency domain. Ad- 
ditionally, within the dipole approximation, both dipoles 
feel the effects of the incident radiation as rebroadcasted 
through electric dipole fields emanating from each other. 
In symbols, that is 



d(i) = aM ■ (Eo + Ep) 
p(i) =5p-(Eo + Em), 



(1) 



where the electric dipole fields Ep = A^'^p^^-' and E^j = 

A^'^d'^^ are expressed in terms of the matrix elements 
A^'^ = [3r^f°' — 6^"']/r^ of the second-rank dipole tensor 
A a — x,y, z). It is assumed here and throughout this 
article that the distance r = |xi — X2I between molecule 
and particle is much less than the wavelength associated 
with Eq. 

In this way, enhanced Raman scattering from the cou- 
pled system of dipoles is achieved by allowing d^^^ and 
p*^^^ in Eq. ^ to couple to and reach self-consistency 
with each other. From the total dipole moment d^^^ -t- 
p(^) = 5tot • Eq, the total system polarizability 



cttot = (1 - OiM ■ A • ap • A) ^ • aM • (1 + A • ap) 
-I- (1 - 5p • A • aM ■ Ay^ • (5p • (1 + A • aM) 



may be directly obtained. The inverse matrix (1 — 5^/ • 
A • ap ■ A)~^ accounts for the classical image effect 
where the molecule's dipole moment repeatedly interacts 
through the particle with its own image dipole, while the 
term A • <5p (recall Ep ^ A ■ ap - Eg) accounts for the 
local-field effect of the particle's dipole electric field upon 



the molecule. The second term accounts for the analo- 
gous image and local-field interactions where the roles of 
molecule and particle are reversed. Both of these effects 
will be generalized to a quantum-mechanical framework 
in the following. 

The enhanced Raman polarizability 



asERS = AQ{d/dQ)at 



(3) 



is related to 5tot by differentiation along the coordinate 
Q with amplitude Q [H, [111 . After some algebra, it is 
found that the enhanced Raman scattering intensity in 
the direction k' with polarization A' from the coupled sys- 
tem of dipoles is composed of the product of two terms: 
a normal Raman scattering intensity from the molecular 
dipole d'^^ in the absence of the particle dipole p'^-*, and 
an enhancement factor stemming from the self-consistent 
coupling of the two dipoles, i.e., 
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^^(A/"' + l)k'A'leA'(k') • SsERS • eA(k)|' 



^Raman(k') , , r , |2 

— enhancement factor , 

4(k) ' ' 



(4) 



where the normal Raman scattering intensity 



Raman (k ) 

/o(k) 



M 



dQ 



^k) 



(5) 

Here, k and k' are the wave vectors of the incident and 
Raman scattered electric fields with associated frequen- 
cies Wk and Wk' and polarizations eA(k) and eA'(k'), 
where A, A' — 1,2. The intensity of the incident field is 
denoted by /o(k) and N'{^,x is the number of Raman- 
scattered photons in the direction k' with polarization 
A'. For the purpose of computing the intensities in Eqs. 
(|3]) and ([5]), we have taken the amplitude of Eg equal 
to i 27r?icjkA/'kA /-^^ca (k) , where is the quantization 
volume of the electric field. This ad hoc choice was taken 
in anticipation of comparison to later results. 

In the limit where the molecular and particle dipole 
moments are aligned with each other and with Eq and 
point along the z-axis, the enhancement factor in Eq. (j4]) 
is greatly simplified and takes the form [sj] 



enhancement factor = 



[1 + (2/r3)5p] 



[l-3p5A/(2/r3)2 
|l + (2/r3)5pf. 



(6) 



(2) where it is assumed that cjk — ^k', otM and ap are 



isotropic, and, in the second line, that ap 



and 



a A/ <C r . This result, which is a straightforward ex- 
tension of the ideas presented in Ref. [3l[ , demonstrates 
the characteristic fourth-power behavior of the SERS en- 
hancement arising whenever a second polarizable body 
is able to react to and act back upon the first: with 
the largest enhancements occurring at a resonance of Sp. 
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Here, it should be pointed out that the SERS enhance- 
ment observed in experiment would scale in the same 
way if both the incident frequency and Raman-scattered 
frequency were the same. Further, it is clear that any 
quantum-mechanical generalization of the SERS mecha- 
nism should reduce in the appropriate limit to the clas- 
sical results reviewed here. 



III. QUANTUM-MECHANICAL APPROACH 
TO SERS 

The quantum many-body problem of SERS, consist- 
ing of a molecular-electronic system adsorbed to a metal 
particle having dimensions between tens and a few hun- 
dred nanometers and containing hundreds of thousands 
of atoms, poses a difficult theoretical and computational 
challenge. Response of such a system to an external 
perturbation presents one example where the associated 
properties are often sufficiently average to permit, to low- 
est order, a first-principles theoretical description [ssj. It 
is precisely this situation that we now study in detail. 

A. Green's functions and many-body Hamiltonian 
of coupled molecule-particle system 

Basic to the theory of response of a zero-temperature 
quantum-mechanical many-body system to an external 
perturbing field is the one-body Green's function ^36|, iBTI . 

MM 

ig{^,t;^X) - «|T{>i'(x,i)*t(x',t')}|0, (7) 

defined as the expectation value of the time-ordered prod- 
uct of Heisenberg field operators 'I' and in the normal- 
ized and interacting many-particle ground state | $^ ) . It 
measures the probability amplitude for an extra particle 
(or hole) to propagate under the action of the full many- 
body Hamiltonian H = h + V through an interacting 
assembly of N particles from the space-time point (x', t') 
to a later (earlier) point (x, t) , there being found in the 
ground state l^^). We note that application of Green's 
function methods to finite, nonuniform systems such as 
molecules [13, SH S!] , atoms [H, |43| , and atomic nuclei 
[ssl . l45j have been extensively studied in the literature; 
see, e.g., Ref. [i^ for a modern account. 



Here we focus on the many-body SERS problem con- 
sisting of a Raman-active molecular system coupled to a 
nearby metallic particle supporting collective oscillation 
of its conduction electrons driven by an external perturb- 
ing field. As we are interested only in the dynamics of 
the molecular electrons in this article, we henceforth spe- 
cialize to the case of fermion statistics where ^ and '^'^ 
are electron field operators satisfying standard anticom- 
mutation relations at equal times, and is the exact 
TV-electron ground state normalized to unity. Effects of 
electronic spin do not alter our proposed methodology in 
any significant way; we therefore omit the spin degree of 
freedom from all equations for simplicity. 

In order to describe the quantum-mechanical interac- 
tion between molecular electrons and a metal particle 
illuminated by an external electric field, it will be neces- 
sary to introduce additional quantum fields for photons 
and for the quantized collective excitations of conduction 
electrons in the metal, known as plasmons. However, 
it is well known that the basic features of the electro- 
magnetic enhancement mechanism of SERS can be ex- 
plained by a classical description of the metal and field; 
see, e.g., Ref. [l6|. In light of this fact, it is our desire 
to build a theory where both the conduction electrons 
in the metal particle and the external field are treated 
classically, yet the molecular-electronic system remains 
quantum mechanical. Henceforth, we assume that the 
quantized electric field Eq is well approximated by the 
classical field Eq {t) — Eq"^' {t) + Eq~ V|) , and we further 
assume a Bogoliubov decomposition [47[ of the quantized 
plasmon field 

il(x, t) = y/n{x,t) + 5$l(x, t) (8) 

in terms of the classical conduction electron density n 
and small-amplitude quantum fluctuations SCl around n. 
The specific form of this classical component will be ex- 
pounded upon in the following. We note that J n{x,t)d^x 
is approximately equal to the total number of conduction 
electrons participating in collective excitation; nonethe- 
less, such excitations will continue to be called plasmons 
in spite of their representation by a classical field. 

In terms of ^' and Cl, the many-body Hamiltonian of 
the interacting molecular electron-plasmon assembly may 
be written in second quantization as 



iJtot ^Vn+ /*^(x,t)/i(x)*(x,t)d3a:-|- /17^(x,<)/t(x)$7(x,<)d^a; + (l/2)/*^(x,i)*^(x',<)y(x,x')#(x',t)#(x,t)d3a;fi^x' 
+ (1/2) / Ot (x, t)Cl\x', t)V{x, x')0(x', t)Cl{x, t)d^xd^x' + J *t (x, i)f2t (x', t)V{x, x')Cl{x', t)#(x, t)d^xd^x', 

(9) 

I 

where a Born-Oppenheimer separation has been made between electron and nuclear coordinates [1^ and all off- 
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diagonal terms in the field operators are omitted. Here, 
Vat = VAr({R}) is the nuclear- nuclear repulsion energy 
which depends upon the set {R} = {Ri,...,RAf} of 
all M nuclear coordinates; note that the nuclear ki- 
netic energy Tjv = TAr({R}) does not appear in Eq. 
^ a& H describes only the molecular electron-plasmon 
physics. Terms two and three, involving /i, represent the 
single-particle Hamiltonians for noninteracting electron 
and plasmon systems, including the external electron- 
nuclear attraction, perturbed by a common externally 
applied classical electric field Eq within the dipole inter- 
action approximation. Introducing molecular and par- 
ticle dipole moments d{t) — J *I'^(x, i)(-ex)^(x, 
and p{t) « /(— ex)n(x, respectively, this perturba- 
tion takes the standard form — d • Eg and — p • Eq, where 
the electric field is evaluated at the molecular-frame ori- 
gin (xi — 0). Terms four and five represent the po- 
tential energy of pairwise repulsive interaction between 
electrons and, separately, plasmons; the last term ex- 
presses the potential energy of interaction between molec- 
ular electrons and plasmons. The two-body potential 
T/(x, x') = e^/|x — x'l appearing above is the instanta- 
neous Coulomb interaction with electronic charge — e. 



We now invoke a second Born-Oppenheimer-like sepa- 
ration between molecular electrons and metallic conduc- 
tion electrons; here we consider only the molecular part 
of the Hamiltonian in Eq. ^ including its interaction 
with the external radiation field and metallic conduction 
electrons. Effects stemming from the noninteracting plas- 
mon Hamiltonian (term three) as well as the plasmon- 
plasmon interaction (term five) are implicitly accounted 
for through the dynamics of the particle's induced elec- 
tronic density, which evolves under the action of 

Hp = W + Jn^x,t)h{x)n{x,t)(fx 

+ (i/2)/r2t (x, t)n^{x\ t)viK, x')f2(x', t)n{x, t)d^xd^x', 

(10) 

where W, defined below in Eq. ([23]) . is the interact- 
ing molecular ground-state potential energy dependent 
upon the underlying classical plasmon density n. The 
associated response of the particle is detailed below in 
Sec. IHIE3I Eliminating these terms leaves the simpli- 
fied molecular Hamiltonian 



H = Vn + J *t (x, t)/i(x)^'(x, t)d^x + J V{^, x')*^(x, i)[(l/2)^'t(x', t)*(x', t) + Cl^^' , i)f2(x', t)]*(x, t)d^xd^x' 
~Vn+ J^-^Cx, t)/i(x)^'(x, t)d^x + /V^(x, x')«'^(x, i)[(l/2)^'t(x', t)«'(x', t) + n(x', i)]#(x, t)d^xd^x', 

(11) 



where the plasmon field operator is expanded to lowest 
order and the products wCt are replaced by the classical 
conduction electron density n. Here, h = ho — (—ex) ■ Eq 
and ho contains the molecular-electronic kinetic energy 
and external electron-nuclear attraction. 

The Heisenberg electron field operator >I'(x,i) ~ 
e:iqy{iHt/h)i'{x.)ex.p{-iHt/h) with *(x) = I]p</'p(x)cp 
may now be expanded onto the basis of electron annihi- 
lation operators 

Cp = (1 - pp)dp + ppbl (12) 

expressed in terms of the particle and hole operators 
dp and bp through canonical transformation. These 
fermionic operators all satisfy standard anticommutation 
relations. The occupation numbers pp, which are eigen- 
values of the exact one-body reduced density 

p(x,x') = «|*t(x)*(x')|0 =5]0;(x)pp,0,(x') 

pq 

(13) 

with matrix elements ppq = \cj,Cq\^o ) ^ take values 
of or 1 depending on the occupation of the pth single- 
particle state. 



In terms of these basic operators and one-body wave 
functions (molecular orbitals) 0p(x) — (x|0p) = (x|p) — 
(x|cj^|vac) with physical vacuum |vac), the electronic 
Hamiltonian reduces to 

H = Vn+ Y,{p\h + U\q)clcg + IY1 {pq\V\rs)clclcsCr, 

pq pqrs 

(14) 

where the sum runs over all N electrons of the molecu- 
lar system, and the explicit time dependence of the ex- 
ternal fields have been omitted for simplicity of nota- 
tion. The molecular dipole moment is expanded as d = 
J2pq{p\ ~ s^\Q)CpCq and included within the one-electron 
Hamiltonian h. Here, U{x.,t) = J V{'x.,x')n{x' ,t)d'^x' is 
the plasmon potential with matrix elements 

{p\U{t)\q) = Jd^x'{p\V\q){^')n{^',t). (15) 

The collective electronic density of the particle 

n(x, i) = no(x) + (5ri(x, i) (16) 

may be decomposed into the sum of a static density 
no and low-amplitude excitations 5n = riint induced 
by some interaction to be specified later. Similarly, 
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U{x, t) = (7o(x) + /7int(x, t) factors into a static potential 
Uq and an induced potential Umt associated with uq and 
TT-int respectively. 

Application of Wick's theorem reduces the above elec- 
tronic Hamiltonian (fH|) to the sum of three parts: a 
scalar constant 

Eo = VN+'^{k\ho + ex-Eo + U\k) + ^^{kl\V\kl^lk), 

k kl 

(17) 

a one-body term quadratic in electron operators 



ground state vector |<&^) stemming from the interacting 
molecular Schrodinger equation 



H\^'^ } = [So + Fo + Fint + HintW^ ) = W^|$o )■ (23) 

The ground-state energy W depends parametrically upon 
the underlying classical plasmon density n. Separation of 
H into an unperturbed component Hq ^ Eq + Fq and a 
perturbation Fiat + Hmt provides a natural ansatz for the 
application of many-body perturbation theory (see, e.g., 
Ref. [H). This wiU be the subject of Sec. HITCl 



P = Y. KpI'^o + ex • Eo + ;7|g) + {pk\V\qk - kq)] N{clcq} 



pqk 



h + ^(p|ex ■ Eo -I- Uiat\q)N{clcq}, 

pq 



(18) 



and a two-body term involving the product of four elec- 
tron operators 



1 



pqrs 



(19) 



In the above, normal ordering is taken with respect to 
particles and holes (i.e., Op and hq) rather than basic 
electron operators Cp. 

The quadratic term F in Eq. (fT8|) is decomposable into 
the sum of the Fock operator Fo — '^pq{p\FQ\q)N {c^pCq} 
and a one-body electron-field and electron-plasmon in- 
teraction 



Fw 



E 

pq 



(p|ex-Eo + C/intk)iV{4cg}. 



(20) 



Choosing the underlying molecular orbitals (pp to di- 
agonalize the Fock matrix {p\FQ\q) = (p|ft.o + Uo\q) + 
'^k{pk\V\qk — kq) = SpqS^, yields the famihar HF orbital 
equation with effective potential Uq 

[/io(x) + C/o(x)]</),(x) 



B. Interaction picture 

We now adopt the interaction picture with respect to 
the uncorrelated noninteracting iV-electron HF ground 
state |$Hp} — ripCpl'^ac) at zero temperature; it sat- 
isfies Fo|<i>j^p) ~ ^SfeSfel'^'Hr) normalized to 
unity. The reference state I^'hf) Fermi vac- 
uum for the electron creation and annihilation operators 
(and particle- hole operators). Henceforth, all expecta- 
tion values will be computed within |<&Hp). The molec- 
ular orbitals (f)p underlying |$Hp) are determined self- 
consistently and satisfy (jp\q) — Spq; it will assumed that 
both |<i>Hp) and the corresponding orbitals and orbital 
energies are known. 

The zero-temperature HF one-body Green's function 
[ssl [soj is defined with respect to the Fermi vacuum as 



$^Ip|r{*(x,<)*t(x',t')}|$^F) 

(/)p(x)zG(",)(i,t>;(x'), 



E 

pq 



(24) 



where iG'^q\t,t') = {^^F\T{cp{t)clit')}\^^p) represents 
the Fock-space matrix elements of iG'^"-' and Cp(t) = 
e:xjp{iFot/h)cpexp{—iFot/h) = exp{—iept/h)cp in the 
interaction-picture representation. The orbitals (f)q are 



solutions of the HF equations (PT|) . In Fourier space, it 
has the Lehmann spectral representation 



+ 



k 



:(fc|F|fc)(x)</),(x) - {k\V\q){x)M^)] = e%{x) G(°)(x,x';c.) ^ ^ 0p(x)G(",) (^)0;(x') 



(21) 



for (fiq and with orbital energy ejj — {q\ho + Uo\q) + 
'^i;{qk\V\qk — kq). Here, and throughout, (pg|l/|rs — 
sr) = {pq\V\rs) — {pq\V\sr) are antisymmetrized matrix 
elements of the two-body potential V. 

Through this normal ordering, the electronic Hamilto- 
nian in Eq. (fTTjl now takes the M0ller-Plesset form [4§| 



pq 



(25) 



H — Eq + Fq + Fi 



int ; 



(22) 



where Fq = Ep^^^i^^p} = Ta^l^l^a - Ek^blk 
is diagonal when expressed in the HF basis. Associated 
with it are the formally exact ground-state energy W and 



where the eigenvalue pp of the HF one-matrix ppq ~ 
(*^HFlCpCg|$Hp) = PpSpq is the occupation number of the 
pth single-particle state and 0"*" is a positive infinitesimal 
needed only to damp and subsequently converge certain 
Fourier integrals; 0+ should be taken to zero at the end 
of all calculations. The retarded component of the time- 
ordered Green's function G*^"^ 



^0^ 



(26) 
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will be of use in the following. In Eq. p4| . Gpq is the ze- 
roth order propagator for an extra electron to propagate 
from the state q to the state p within the A^-particle non- 
interacting background I'&Hp}. In this approximation, the 
electron does not interact directly with any other particle 
nor is it scattered out of the single-particle state p (due 
to the delta function 5pq in Eq. (f25|) ). Interaction effects 
will now be systematically built in through the machin- 
ery of time-dependent many-body perturbation theory 

C. Many-body perturbation theory 

Until this point we have developed an essentially exact 
many-body theory of SERS including the coupling of a 
general molecule to a plasnion-supporting metallic sys- 
tem under the influence of an external perturbing elec- 
tric field. Building from the noninteracting HF refer- 
ence state of the molecule, it will now be demonstrated 
how to systematically incorporate these interaction ef- 
fects to arbitrary order through the machinery of many- 
body perturbation theory. As a first step in this direc- 
tion, we choose to consider only those interactions origi- 
nating from the external electric field and from collective 
excitations of conduction electrons in the metal particle 
induced by both the external field and by interactions 
with the molecule, and truncate the level of molecular- 
electronic correlation to HF mean-field theory; there- 
fore, we suppress the two-body electron-electron interac- 
tion _ffi,it in Eq. (|19p and focus upon the perturbations 
stemming from -Flnt- It is, however, important to note 
that, due to the generality of our approach, electronic- 
correlation effects (from i?int) may be included by using 
the same perturbation techniques described below, or, al- 
ternatively, may be added in the spirit of a Kohn-Sham 
density- functional theory approach . Within this con- 
text we note that the chemical mechanism of SERS may 
be explored by including a subset of the metallic con- 
duction electrons in addition to the molecular electrons 
within the density- functional theory; see Appendix \K\ for 
details. 



1. Interaction of molecular electrons with a metal particle 

We focus first on the perturbations of the molecular- 
electronic system induced by the presence of a metal par- 
ticle nearby. From Eq. pB]) . the collective electronic 



density of the particle 

n(x, i) = no(x) + nji/(x, i) (27) 

may be decomposed as the sum of a static density tlq 
and low-amplitude excitations 5n = hm induced by 
the molecule itself. This will be shown to underlie the 
quantum-mechanical analog of the image effect discussed 
in Ref. [sT'l . Associated with these densities are the static 
and induced dipole moments po and Pm of the particle. 
With this decomposition, the interaction between molec- 
ular electrons and the induced density um in the particle 
can be written to lowest order as 

F^fit) - F/^(t) = {p\Jd^x'VMnM{^',t)\q) 

^-(p|-ex|g).^^[P^(^) J"^"^^'^ 

(28) 

following multipole expansion of the potential V, where 
xi — X2 = rf and where the pa/(0 = /(— ex)nA/(x, t)d^x 
is the dipole moment induced in the particle by the 
molecule. Effects of retardation are neglected. In prin- 
ciple there is no reason to additionally impose a dipole 
interaction approximation between molecule and parti- 
cle; however, since we are interested here in developing a 
quantum-mechanical generalization of the classical dipole 
model presented in Sec. |TT1 we choose to do so. Note 
that the interaction between molecular electrons and the 
static part of the density uq is already accounted for in 
the reference system; see, e.g., Eq. pT|) . The above ex- 
pression represents the interaction energy between the 
molecular dipole moment d and a molecule-induced ex- 
citation of conduction electrons in the metal particle. It 
can be rewritten in more compact form as 

Fint(t) ^ Mt) = 5](p|C/p(i)|g)Af{ct(t)c,(i)} 

pq (29) 
« -d(t) • A-pM(t), 

where, as before, A^°' = [Sr^f"^ — d^'^]/r^ are the matrix 
elements of the second-rank dipole tensor A. 

Dyson has provided an algorithm for computing the 
exact interacting iV-body Green's function by perturba- 
tive expansion from the noninteracting reference state 
of the interaction picture f39l|, where the effects of in- 
teraction are encapsulated within the irreducible self en- 
ergy ftE* [sij l . In the particular case of the perturbation 
Fp = — d ■ A • Pm, Dyson's expansion yields 
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{-i/nr 



Jdh--- dtn{'i>^F\T{[~d{h) ■ A • pMih) d(t„) • A • PM(t»)]Cp(t)ct(t')}|*HF)c 



n=0 



^5]/dVdMG(";(t,t;)(r|ex«|t)A«'^(ex?)jnp(xi,ii;x;,t;)(exf)A^^Gl^^ 



r 



(30) 



where (■ ■ ■ )c indicates that only Unked or connected di- 
agrams contribute in the expansion [52| . 

The classical polarization propagator of the molecule- 
induced density fluctuations um = n — hq in the particle 
can be identified in the second-order term above as 

znp(x,t;x',t') =nAf(x,t)nM(x',t') (31) 

in analogy to the quantum-mechanical noninteracting HF 
polarization propagator [H, [s^ defined in Eq. (|B1|) in 
Appendix [b1 From the retarded component of Hp, the 
molecule- induced classical linear polarizability of the par- 
ticle can be written as 



-iha^p{t,t') ^ J(fx(fx'{-ex^)iU^{x,t;y:',t'){-ex"') 
4( 

(32), 



and has a mathematical structure similar to the noninter- 
acting HF molecular linear polarizability (jB3[) defined in 
Appendix [BJ Following omission of the first-order per- 
turbative correction, which is of no importance in the 
following (and can be renormalized away), the retarded 
component of Eq. (15^ becomes 



E. 

rs 



J dhdt[[Gfo-^]pr{t, ii) [— ^(r|ex|t) • A • ihapih,t[) ■ A ■ [Gfo)]t„(t'i, ti)(u|ex'|s) 



[G(o)]s9(ii, i') + 



(33) 



with irreducible self energy 
;i[S*U(ii,tl) = ^^(p|ex|r).A 

rs 

■ ihap{h,t[) ■ A • [Gfo)],,(t'i,<i)(s|ex'|(?) + • • ■ . 

(34) 

Here, the first term on the right hand side already con- 
tains the desired effects of polarization. We, therefore, 
truncate the above perturbation series for /lE* at sec- 
ond order in Fp and define the spectral representation of 
the second-order irreducible self energy by 

;.[E*](f)(c.)=^(p|-ex|r).A 

rs 

duj' — 
■ J — 5p(w')[Gfo)]rs('^ + uj')-A- {s\ - ex'\q). 

(35) 



An approximate yet infinite-order nonperturbative rep- 
resentation of the exact interacting one-body retarded 
Green's function may now be constructed by resumming 
the Dyson series 

= [Gfo)]p,(.;) +5][Gfo)]p.(.;)p*](f^2)(^)[g^i,,(c.) 

rs 

(36) 

stemming from the second-order perturbative approxi- 
mation S* K, [5^. Such an approximation cor- 
responds to the physical scenario where the molecule is 
able to repeatedly excite density fluctuations in and po- 
larize the particle and these excitations act back upon the 
molecule to inflnite order. Said differently, the quantum 
analog of the classical image effect discussed in Ref. 
is included by resumming the Dyson series based upon 
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the above second-order truncation of the self energy. For 
simphcity of notation we drop the label 2 in the follow- 



ing so that [S* 



p — 



Note that [S*]^ contains 



only one of several irreducible terms occurring at sec- 
ond order in perturbation theory (39j : however, among 
all others, this is the only contribution that includes the 
desired polarization effects. 

In similar spirit, the approximate time-ordered one- 
body Green's function may be compactly expressed 
by the Dyson expansion 



(37) 



In analogy to Eq. (|36p , solving this recursive equation for 
resums the infinite class of diagrams spanned by the 
second-order perturbative truncation of the time-ordered 



self energy E* s 
inverse solution 



P*](p,2) = J.* Equation ^ has the 



[GpXqiu;) = [Gl,]U{iu) - [J:*U{u;), (38) 



which may be written in matrix form as 



(o)J' 



■'p> 



-'Pv 



(0)J 



(39) 



where we have used the fact that the one-body HF 



(o)J» 



Green's function ((25|) is diagonal in its indices; [G 
Sp, and [G^^]o — H*^ represent diagonal particle- 
particle and hole-hole matrices, while S*,^ and S*v rep- 
resent off-diagonal particle-hole and hole-particle block 
matrix contributions to Qp^- Formal inversion may be 
expressed in terms of minors as 



= [detgpi]-i(-)P+«minor[api],p. (40) 



Alternatively, we make use of the block-matrix inverse 



([G 



([G 



(0)] 



(o)J 



^Po) ^^Pv([G 



(O)J' 



]. - )- 



([G 



(0)J 



([G(o)]o - Sp^j) ^ 



■'Po) 



(41) 



together with the weak-coupling approximation where which are here derived in the weak-coupling limit. 
[E*]p,,[S*]gp <C [E^]pp, [E^],,. Here, the infinitesi- 



mals ±iO+ are omitted from G'-'^-' as the self energy 
(I35p is complex valued and has an imaginary compo- 
nent that changes sign with respect to the chemical po- 
tential according to ImEp(w) > when u < ii/h and 

ImE*(w) < when uj > n/h. Like G^°\ in Eq. gl|) 
accounts for the propagation of an extra electron from the 
single-particle state q to p] however, unlike the noninter- 
acting Green's function where p = q, Gpq also permits 
electrons to scatter out of q and into the single-particle 
state p ^ q. This latter process is described through the 
off-diagonal components 



This interacting molecular-electronic one-body Green's 
function has the space- and frequency-dependent 
form 



g^(x,x';c^)= Vxp(x)ap^,(a;)x;(x') 



pq 



(43) 



-C.~£0/ft-[E*],,(c.) 
Sab 



Sab 



w-eg/n- [E*]f,b(cj) 



with interacting Dyson orbitals Xg(x) = ^^0r(x)C/rg 
satisfying the nonlinear integro-differential equation 



^p]bt(^) 



-^c.-£0/;i-[E*],,(c.) ^ L,-e]lh-{Y.%M 

2) 



[/lo(x) + C/o(x)]x,(x) +^[(fc|y|fc)(x)x,(x) - (fc|F|g)lx^Xfe(x)] +/dV;iE*(x,x^^)x,(xO = + ^^^^ 

fe 



expressed in terms of the nonlocal frequency-dependent (energy-dependent) effective potential Ep(x, x';cl;) 
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EpgXp(x)[Sp]pq(w)x*(x')- It is derived by applying the 
operator 

iHF(x, x'; uj)^[nuj- {ho{x) + C/o(x)}](5(x - x') 

- J2mV\k){^)S{^ - x') - x:.(x')V^(x,x')xfc(x)], 

k 

(45) 

which is defined in terms of Xk, to in Eq. ^ and 
then projecting against Xq', -^hf can equivalently be writ- 
ten in terms of the HF orbitals (j)k where it satisfies 
/ (i3yLHF(x,y;tj)G("ny;x';w) = M(x - x'). Equation 
incorporates the image interaction associated with 
the molecule-induced excitations of the particle into the 
molecular HF mean-field equations pT|) through the ir- 
reducible self energy ftS* in Eq. ([55]) . These interacting 
orbitals and orbital energies, which are solutions of Eq. 
(|33]), must reach self consistency with h['S*]pq defined in 
Eq. ([58|) as the self energy both determines and is de- 
termined by the new interacting orbitals. In this way, 
consistent solutions of the interacting molecule-particle 
system may be obtained. 

2. Interaction of coupled molecule-particle system with an 
external electric field 

Now that we have derived an expression for the molec- 
ular one-body Green's function interacting with its self- 
induced electronic density fluctuations in a nearby metal 
particle {i.e., its image), we are ready to build in the 
perturbing effects of an external classical electric field 

Eo(t) = E[,+^e"*'^''* -I- E[,~^e*'^''* (46) 

with frequency lu^ upon the combined and coupled sys- 
tem. The field interacts directly with the molecule and, 
additionally, indirectly by inducing small- amplitude col- 
lective excitations of conduction electrons in the particle 
described by n^. This latter interaction is described in 
Ref. [Slj as the local-field effect. Following Eq. (|16p . the 
particle's electronic density can be decomposed as 

n(x, t) — no(x) + n£;(x, t) 

= no(x) -|-n^+'(x,t) -|-n^"\x,t) 

in the linear response limit of the external field, where the 
labels (-I-), (— ) refer to excitations set up by the incoming 
and outgoing components Eq^'' and Eq ' respectively of 
Eq. Despite our classical treatment of the external elec- 
tric field, we continue to speak of photons and will, when 
justifiable and appropriately pointed out, need to make 
an ad hoc adjustment of photon occupation number as a 
result of our classical-field ansatz. Of course, we could 



have treated the external field quantum mechanically 
and introduced the appropriate photon Green's functions 
needed to carry its dynamics. However, it is our aim to 
keep this presentation as clear as possible and elucidate 
only the lowest-order processes governing SERS, which, 
save a small error in photon occupation number, is de- 
scribable without resorting to field quantization. 

Multipole expansion of the electron-field and electron- 
plasmon interaction Hamiltonian i^int, given in Eq. (|20p . 
results in the following expression 

F^fit) = -{p\ - ex|g) ■ Eo{t) + (p|/d^Tni;(x,i)F(x,x')|g) 

^-(p|-ex|,).{Eo(t) + ^^[P^(^)^;]-P^(^^ } 

(48) 

for its matrix elements at dipole order. The second term 
in brackets is the classical electric dipole field Ep(i) = 
^ ■ P-e(0 — ^ ' /(— ex)n£;(x, t)(i^x of the metal parti- 
cle's dipole plasmon as induced by Eg. Like Eq. it 
contributes the effective plasmon potential 

C/i„t(x,t) UE{^,t) w -(-ex) ■ A • pE{t) 

= -(-ex).Ep(t), ^ ^ 

where, like ue and pe, Ep is an abbreviation for 
E^-* -t- Ep . Since we choose to specialize to the dipole 
approximation for simplicity, no other multipole contri- 
butions will be considered. Therefore, perturbations to 
the molecular-electronic system from Eg arise directly 
and indirectly through 

Fintit) ^ pEit) - -d(i) • [Eo(t) + Ep(t)] 

= -J2{P\ - e^\q)N{Cl{t)Cq{t)} ■ [Eo(t) + Ep(i)], 

pq 

(50) 

where the new interacting molecular electron operators 
C'p(i) = J2q UpqCq{t) are related to the old noninteracting 
operators Cp{t) by the same unitary transformation U 
that affected the molecular orbitals. 

As before, Dyson's expansion provides a systematic 
way to build in the interaction effects of Fe into the 
molecular one-body Green's function. In contrast to 
the previous application of the Dyson expansion where 
the interaction effects of Fp were included on top 
of the noninteracting HF one-body Green's function 
G'^°\ here we build the interaction effects of Fe upon 
the interacting one-body Green's function Qpq(t,t') = 

($^p|r{Cp(i)C't(i')}l$^p) that describes the coupling 
between molecular electrons and their image as mediated 
by the conduction electrons of a metallic particle. We see 
that 
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E 

n=0 



J dh ■ ■■dtn{'S>^p\T{-d{h) ■ [Eo(ti) + Ep(ii)] • ■ • - d(i„) • [Eo(i„) + Ep(i„)]CpWCi(t')}|*HF)c 



+ ^^/diigp^,(i,ii)(r|ex|s)ef^(ii,t') • [Eo(ti) + Ep(ti) 



fi2 



rstu 

+ Ej,{h)ES{h) + E%{h)E^p{t2)] + ■■■ , 



(51) 



r 



where, as was discussed previously in the context of Hp, 
the time-ordering affects the classical fields through 



i?o^(ii)Kfe) = r{£;«(ii)i?ofe)} 

- B{ti - h)El{h)E^{h) + e{t2 - h)E^{t2)El{h). 



(52), 



Similar expressions can be written for the remaining 
three terms in curly brackets above. Truncating the per- 
turbation series in Eq. ([5T|) at second order in Fe results 
in the second-order perturbative approximation to the 
interacting one-body Green's function 



rs 

rs tu 

+ El{h)E-p(t2) + Eiih)ES{h) + EUti)Ef,{t2)}] gf,{t2,t'), 

I 



(53) 



where we have omitted the first-order perturbative cor- 
rection as it describes only the stimulated absorption and 
emission of Eq and E p which are of no importance in the 
following (and can be renormalized away). 

Indeed it is not until second order in the external field 
perturbation that Raman scattering from the molecular 
system can be described. Other competing processes oc- 
cur at second order as well, such as those of two photon- 
absorption or two-photon emission. However, we are not 
interested in describing these events and, consequently, 
prune away all terms in Eq. (j53p not related to Ra- 
man scattering. Due to the presence of both Eq and Ep, 
several types of Raman processes are present in the re- 
maining expression. Note that the electric field Ep stems 
from collective excitation of conduction electrons in the 
metal particle induced either by the external field (i.e., 
the local-field effect) or by the Raman-scattered field of 
the molecule. In either case, it may be written in Fourier 
space as 



(±) 



{A • ap{uj) ■ E[)^^}27r5(a' ± cjk,k' 
E[,p^(w)27r(5(w±cjk,k') 



(54) 



in terms of the linear polarizability defined in Eq. 
where the wave vector k is associated with the incoming 
field labeled by {+) while the wave vector k' is associ- 
ated with the outgoing field labeled by (— ); they satisfy 

[e(+V = e^^ 

Together with this expression for E^'', the interacting 
one-body Green's function in Eq. ([55)1 enjoys the Fourier 
decomposition 

g'X^'H^) = qU^) + E Q^A^)^t}-^\^)QU^) (55) 

rs 

with second-order irreducible self energy for Raman scat- 
tering 

fiS*(^^2)(..) « \ Y.^\ - ex^\r)g^M{s\ - ex"^\q) 



X < 



EQpli-^'k')E^l' + E^p'^{uji,)E^^ 

^0? ^^o^l('^k) + E^^'^^ e1^pI{-uji,,) 

_ +^(p)(-a>kO^^Jl(^k) +^^S(^k)^^pU-'^k'), 

(56) 
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FIG. 2: (Color online) Enhanced Raman-scattering Feynman diagrams occurring at second order in perturbation theory in 
the interaction Fe = —d ■ [Eq + Ep]. Diagram (a.) is analogous to normal Raman scattering of the field Eq'^^ Diagrams (b.) 
and (c.) represent a particle-mediated Raman-scattering process where either the Raman scattered field [diagram (b.)] or the 
incident field [diagram (c.)] is scattered by the particle. Diagram (d.) represents the process where both the incident and 
Raman-scattered fields are mediated by the particle. This particle-mediated field E^p' = A ■ ap ■ Eq"^' is defined in terms of 
the particle's polarizability ap (|32p . While retained in the formalism [see, e.g., Eq. (|56|l ]. crossed diagrams are not drawn for 
simplicity in presentation. However, all two-photon absorption and two-photon emission processes are omitted here and in the 
theory as are all first-order perturbative contributions. Note that all molecular-electronic propagators are interacting one-body 
Green's functions defined in Eq. (|4ip . 



where the energies of the intermediate electronically ex- 
cited molecular states r and s both include the addi- 
tional energy HiOk (— fi^k') of a single absorbed (emit- 
ted) photon. By replacing by the first pair 
of terms in brackets represents ordinary Raman scat- 
tering from the molecule in the absence of the parti- 
cle. Terms three and four are associated with the mixed 
event in which the incident field directly interacts with 
the molecule while the molecular Raman-scattered field 
scatters off of the particle before detection; terms five 
and six represent the opposite time ordering of terms 
three and four. The last pair of terms are associated 
with the scattering event where the incident field is first 
scattered by the particle. This enhanced field interacts 
with the molecule, which subsequently Raman scatters 
the radiation back to the particle. In the final step, the 
particle rebroadcasts the molecular Raman field to the 
detector. These processes are all summarized in Fig. [2] 
While they are present in the formalism, crossed events 
where a photon is first scattered and, later, a second pho- 
ton is absorbed are omitted from the figure. As we are 
interested only in Raman scattering, two-photon absorp- 
tion/emission processes, which would involve terms like 

^ot^^ot^ or E^'^^E^Q^^ are not considered. Note that 
additional diagrams would be present had we chosen to 
quantize the electric field. Henceforth, for simplicity of 
notation we drop the label 2 so that S* = 2V 



3. Random-phase approximation 

Until this point we have assumed that the collective 
excitations of conduction electrons in the metal parti- 
cle can be represented entirely by their classical den- 



sity n. When a molecular-electronic system is brought 
into the vicinity of a particle it induces density fiuctu- 
ations in the metal that are, to first order, describable 
by — n — hq. For the purpose of demonstrating cer- 
tain properties of our formalism, we here invoke the fur- 
ther approximation that the conduction electrons of the 
spatially anisotropic metal particle are well-described as 
a homogeneous electron gas (or electron plasma) in the 
high-density limit. This approximation becomes appro- 
priate when the spatial dimensions of the particle are 
larger than the mean-free path of its conduction elec- 
trons. For Au and Ag, which are typical SERS sub- 
strates, the mean-free path of conduction electrons is ap- 
proximately 40 and 50 nm respectively. This justifies the 
replacement of the particle's polarization propa gator by 
the random-phase approximation (RPA) result |54l . l5a | , 
i.e., Hp Hrpa — np(l — FpIIp)"^ = IIpK"^, where 
K = 1 — (ilo/oj)'^ is the generalized dielectric function of 
the particle in the RPA. Similarly, the polarizability of 
the metal particle takes on the RPA form 

5^ (c.) ^ al^p^ioj) = -^S^", (57) 



where fig — ^ A.Tie'^np jra is the bulk plasma frequency 
and np = kp/Sn'^ is the density of the free electron gas 
with Fermi wave vector kp = y^2mep/h?. It is, of course, 
a severe approximation to assume that the polarization 
propagator of an arbitrarily sized and shaped particle will 
have the RPA form. (Extension to a damped Drude or 
Lorentz oscillator model with several resonant frequen- 
cies would be straightforward.) Rather the dynamics of 
the metal particle's conduction electrons and their in- 
duced dipole moments pm should be solved for explicitly. 
While this ultimately is our desire and is the subject of 
Sec. nil E 31 below, such a task requires tremendous nu- 
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merical effort and, at this stage, it is only our intent to 
lay out the basic working equations of our model and 
to highlight some of its general results and salient fea- 
tures. Choosing, for this purpose, to temporarily make 
a detour and invoke the RPA provides a physically rea- 
sonable analytical model of the particle's response that 
is sufficiently rich to allow us to do so without having to 



explicitly compute the electronic dynamics of the metal 
particle. 

The self energy ftS*, accounting for the particle- 
mediated interaction of molecular electrons with their 
own image, is hereafter approximated in the RPA (j57p 
as 



^{p\ -ex|s) - A- {s\ -ex'lg) 



(58) 



eO/h-no-uj + iO+ eyn + no-uj -iO+. 



+ 




FIG. 3: (Color online) Second-order (second Born approxi- 
mation beyond HF) one-body Green's function account- 
ing for the interaction of molecular electrons with collective 
excitations of conduction electrons in a nearby metal par- 
ticle described by the RPA polarization propagator IIrpa. 
Here, the molecular-electronic propagators are one-body HF 
Green's functions G'*^'. 



Here, the near-idempotency of A'''' 



been exploited to simplify the expression. The diagram- 
matic representation of Q^^ is displayed in Fig. [3l Note 
that the sum runs over the set of all single-particle states 
which depend parametrically upon the M nuclear coor- 
dinates. 

The complex-valued and frequency-dependent irre- 
ducible self energy hYi* may be decomposed into real 
and imaginary components as 



[Y.*f{huj + = A^ihio) - {i/2)rP{hw). (59) 



By appeaHng to the identity (x ± iO~^) ^ = P{x ^) =F 
T:iS{x) with principle value P, its real diagonal matrix 



elements 

A;(£°)^Re[SX(£p) 



-P^(p| — ex|s) • A ■ (s| — ex'lp) 



(60) 



account for the shifting of the pth orbital energy while 
its imaginary diagonal matrix elements 



^(i/2)r,^p(4) 

2 



2mQor 



{p\ — ex|s} • A • (s\ — ex'lp) 



X [(1 - pO)<5(£^ - £° + hn,) - pO<5(£^ -el- nno)] 

(61) 

account for the rate of spontaneous emission of a plas- 
mon with energy hD,Q from the electronically excited state 
p (first term) or the rate of spontaneous absorption of 
a plasmon with energy hflo into the state p (second 
term), both inducing molecular-electronic transitions to 
the state s, where we have have assumed that and P^ 
vary so slowly with energy that we may choose huj — 
both effects are due to the interaction between molecu- 
lar electrons and the induced plasma density in the par- 
ticle and are, here, rigorously included from first prin- 
ciples. From the point of view of the molecule, these 
interactions underlie a state-by-state broadening of the 
molecule's electronically excited states. It is precisely 
this interaction physics that is not explicitly treated in 
Rcf. 22, 24, 25], but is, rather, implicitly encapsulated 
within a common phenomenological damping factor for 
all electronically-excited states. A generic consequence 
occurring whenever P^ 7^ is that the effective Hamilto- 
nian of the coupled molecule-particle system is no longer 
Hermitian. For completeness we note that both and 
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are related to each other by Hilbert transformation 

This approximation of the metal particle's electrons as 
a homogeneous electron gas supporting collective excita- 
tion at the bulk plasma frequency is applied in the fol- 
lowing to analytically demonstrate an enhanced Raman 
scattering from the coupled molecule-particle system in 
analogy to the classical result [3l| reviewed in Sec. HIl 




D. Scattering T-Matrix 

Transition amplitudes between initial and final eigen- 
states of an arbitrary reference Hamiltonian underlie 
the computation of many different observable quantities; 
such amplitudes are directly related to the scattering T- 
matrix which is the subject of this section. Recalling that 
time-dependent Green's functions are propagators in the 
sense that they describe the propagation of particles in 
time through an interacting many-particle assembly, it is 
not surprising that a connection exists between the exact 
one-body Green's function and the 5- matrix of scattering 
theory; see, e.g. Ref. Specifically, their relationship 
for t > is given by 



5f- - lim p*(^p+^'')*/2'*f —G^(z) 



= Sn 



(62) 



in the interaction representation, with one-body inter- 
action V, where the effects of scattering from an initial 
many-body state with underlying one-body states labeled 
by (7 to a final many-body state with underlying one-body 
states labeled by p are encapsulated in the one-body T- 
matrix elements [53| 



{p\V\q)+Y,{p\V\r)g^Ae,)mq)- (63) 



Here the retarded Green's function Q^{z) is analytically 
continued away from the real axis and into the complex 
z-plane where new features such as complex poles (also 
called resonances) and complex thresholds may be re- 
vealed on higher or lower Riemann sheets. 

As we will be concerned with perturbations stemming 
from an external electric field Eq, we make the dipole- 
interaction approximation and take V = — d • Eq for 
the purposes of the present discussion. The contour C+ , 
which is displayed in Fig. 01 is rerouted to avoid the 
branch point resulting from the (real) threshold where 
an electronic continuum channel opens due to the ab- 
sorption/emission of a photon from/to the field Eq. An 
associated branch cut connects this branch point to its 
terminal branch point chosen ai z — (-l-oo, 0). As a result 
of this particular route for C+, the contour integration 
moves onto the second Riemann sheet on the right-hand 
side (shown in red in Fig. [¥]), where it encloses resonances 
at the points zi, Z2, and z^ (assuming, for the purpose 
of demonstration, that only three exist) . What were real 



FIG. 4: (Color online) Complex contour associated with the 
S- matrix integral in Eq. (|62|l . The contour wraps around the 
branch point, where a molecular-electronic continuum channel 
opens due to interaction with the external field Eq, and passes 
onto the second Riemann sheet on the right-hand side; the 
associated part of C+ is indicated with a dashed red line. 
There it encloses complex poles (resonances) at the points 
zi, 22, and 23. A branch cut (blue) extends from the first 
threshold to z = (-|-oo,0). 



eigenenergies of the unperturbed Hamiltonian now be- 
come complex eigenenergies (Imzj 7^ 0, j = 1,2,3) of 
the interacting system. The particular locations of these 
complex poles of the exact one-body Green's function are 
due to the analytic continuation of from the upper- 
half plane onto the second Riemann sheet in lower-half 
plane. 



1. Normal Raman scattering from a noninteracting 
molecular system 

We are now in a position to compute Raman transition 
amplitudes between states of the interacting molecule- 
particle system induced by the perturbation V = Fe = 

d ■ [Eq + Ep]. However, before computing the asso- 
ciated interacting scattering T-matrix, we first make a 
detour and consider the case of normal Raman scatter- 
ing from an isolated Raman-active molecule using the 
many-body Green's function formalism. A more thor- 
ough theoretical development of Raman scattering that 
does not involve Green's functions can be found in Ref. 
[5?! , while an advanced review covering linear and nonlin- 
ear optical processes from a Green's function perspective 
can be found in Ref. ^58*1. Here, noninteracting molec- 
ular electrons are described at the level of HF mean- 
field theory by the one-body HF Green 's function G'-°^ 
defined previously in Eqs. (|24p and molecular- 
nuclear degrees of freedom J = 1 , . . . , M, which underlie 
all electronic states within the Born-Oppenheimer ap- 
proximation, are represented by the one-body vibrational 
states \ vj). In the harmonic approximation, the total AI- 
body molecular-vibrational wave function is equal to the 
unsymmetrized (Hartree) product nj=i(Qj|^j) of 
body vibrational wave functions (Qj|i^j} for each degree 
of freedom. Molecular rotations are not resolved in our 
presentation. Perturbed by the external field Eg, the 
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scattering T-matrix of the molecular system may be ap- turbation V = — d • Eg as 
proximated at second order in the dipole-interaction per- 



/ / |^HF(2)| y 



{p\ex^\r)5rs{s\ex"^\q) {p\ex'^\r)5rs{s\ex'^\q) 



(64) 



which is expressed in terms of the retarded HF Green's 
function (|26p . Here, in addition to the sum over 
the intermediate electronic states r and s, there should 
be a sum over the intermediate vibrational states of the 
molecule; however, for simplicity in presentation, we omit 
the vibrational energy differences in the denominator 
(in comparison to the electronic energy differences) here 
and in the following and appeal to the closure relation 
1 = in the numerator [STj . It is also im- 

portant to note that, due to Kronecker delta 5rs in the 
numerator, the intermediate electronic states must be the 
same and, further, must label either particle-particle or 
hole-hole states; no particle-hole or hole-particle interme- 
diate states contribute to Raman scattering. This point 
will be important in deriving an expression for the SERS 
intensity in Sec. IIII D 21 below. 

The one-body states which underlie the initial and final 
molecular states for normal Raman scattering are 




FIG. 5: (Color online) Diagrammatic representation of the 
processes occurring in normal Raman scattering from a 
molecular target with electronic transition polarizability 5jf . 
Molecular vibrational states (red) are labeled by vj and u'j 
and the incident and Raman-scattered electric field (green) 
by Eq^' and Eq~'. Both are implicitly accounted for in our 
semiclassical approach through the state labels. 



\<1)Wj) = \4>q]{vj\M],xML'X'})Wj) . 

|p)|i.,';)^|0p;K;(AA-l)kA,(AA' + l)k'v})|i'S), ^ ^ 



where the incident and Raman-scattered fields are implic- 
itly labeled in the molecular-electronic states only to mo- 
tivate proper field normalization; here A/" (TV— 1) photons 
are initially (finally) in the state |k; A) with wave vector 
k, polarization A, and energy hw\^, and TV' (TV' -I- 1) pho- 
tons are initially (finally) in the state |k';A') with wave 
vector k', polarization A', and energy Swk'. As previously 
discussed, within the Born-Oppenheimer approximation, 
the molecular-electronic and nuclear coordinates are sep- 
arated as \q] {vj})\vj) and |p; Wj})\v'j), where vj and v'j 
label the initial and final vibrational quanta associated 
with the particular normal-mode coordinate Qj. 

Recognizing the inverse of the retarded HF polariza- 
tion propagator H^-j [Eq. (|B2p ] in the denominator of 

Eq. (l64l) and recalling the connection between H|^^ and 



the linear polarizability defined in Eq. (|B4p . we find that 



/ / |^HF(2)| > 



-E 



(-) 



■0 ■{v'J\a%[h^^,,~huJ^,)\uJ)■^'^^ 



(+) 



^ ^.^^/ 27rr^Jk'(AA'-H)k'v ^.^ /27rfiiJkTVkA 



iT3 



L3 



X 47^(k') . {v'j\a'^{{hi,^,-tKo^,)\vj) ■ 4+)(k), 



(66) 



where, as in Sec. |TT1 the incident electric field ampli- 
tude is i^j2'Khu}]i^N-\iix/ L-^e^^^ (k) as is consistent with the 



field occupation numbers in Eq. (j65[) . For completeness 
we point out that had Eq been properly treated as a 
quantum-mechanical field, the number of scattered pho- 
tons in Eq. (p^ would have rigorously been (TV' + l)k'A' ; 
see, e.g., Ref. [s^. A diagrammatic representation of 
the electric field interaction and nuclear vibrational pro- 
cesses occurring in normal Raman scattering is displayed 
in Fig. m 
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Fermi's golden rule of time-dependent perturbation 
theory ^5Q] dictates that the rate of transition w between 
states / and i is related to the transition amplitude by 



w 



HF(2) 



27r 



(67) 



for the particular case of normal Raman scattering from 
a noninteracting molecular system, where p(fiwk') — 
[L /2ncY' {lo^, / K)dQ,\^i is the density of states of the emit- 
ted electric field propagating in the k'-direction. It is now 
straightforward to compute the Raman-scattering inten- 
sity in the direction k' from a single molecular scatterer 
to be 



-^Raman(k ) 



-(AA' + lw|^4'^(k') 



• KKQj - Qo) • VQ,,5;i(fia;k, -nw^')W.j) ■ 4+^(k) 

(68) 

where /o(k) is the intensity of the incident field in the di- 
rection k with polarization A, and where only the linear 
term in the Taylor expansion of the electronic polariz- 
ability 

OLM [{Q}] = 5m[Qo]+^(Qj-Qo)- 

(69) 

around the molecular equilibrium geometry Qo was re- 
tained. Here the sum runs over all M normal-mode co- 
ordinates {Q} of the molecule and it is simple to show 



that {v'j\{CIj — Qo)|'^j) is nonzero whenever v'j = vj± 1. 
Elastic Rayleigh scattering is described by the first term, 
while inelastic Raman scattering, at lowest order, oc- 
curs through the second. Raman-scattering overtones be- 
yond the fundamental depend upon higher-order terms. 
This expression is the quantum-mechanical analog of the 
classical expression derived in Eq. ^ for a polarizable 
molecule interacting with the external field Eq. 



2. Enhanced Raman scattering from an interacting system 



Having briefly reviewed the theory of normal Raman 
scattering from a noninteracting molecular-electronic 
system within the many-body Green's function formal- 
ism, we now turn to the case where the molecule is itself 
interacting with a nearby metal particle, including both 
image and local-field effects. The one-body Green's func- 
tion developed in Sec. IHI C II was designed specif- 
ically to incorporate this physics; it includes the self- 
induced polarization effects of a nearby classical metallic 
particle to infinite order in perturbation theory (i.e., the 
image effects) and, for the purpose of demonstration only, 
assumes that the conduction electrons of the metal are 
well-described by the RPA. Comparison of the expres- 
sion (|63p for the scattering T-matrix with the second- 
order perturbative approximation to the irreducible self 
energy /tS* displayed in Eq. ([56]) shows that 



vs:XY 



E*, ^ • (j>\e^\r)[gXs{el + hi,^,){s\e^\q) ■ E^+^ -t- eV+' • {p\e^\r)[gX,{e'l - n^k')(s|ex'|g) • e|^ ^ 



(+) 



(70) 



which is defined in terms of the retarded interacting one- 
body Green's function Qp^ and where the energy of the 
incident field is huj\^ (60j . The underlying one-body states 
associated with enhanced Raman scattering are similar 
to those previously defined in Eq. (|65p for normal Raman 
scattering, i.e., 

= \xq\{v.rM\.x,K'x'})Wj) 
b)K)^lxp;W;(A^-i)kA,(AA' + i)k'A'})l^:7)- ^ ' 

As previous, they label both electronic and nuclear vi- 
brational states of the molecule. However, here, the one- 
body electronic states are described by the interacting 
Dyson orbitals Xq, which are solutions of the Dyson equa- 
tion ([33]), rather than by the noninteracting HE orbitals 
(pq. The labels X,Y € {0, OP} refer to either the exter- 



nal field Eq^-* or to the field of the particle Epp''. From 
the two terms in Eq. ([70)1 there are eight possible ways 
to arrange these two fields between the two states X, Y : 
four from the first term (uncrossed interactions) and four 
from the second term (crossed interactions). All eight 
terms are included in the formalism; they enumerate all 
possible time orderings between uncrossed and crossed 
interactions. In this sense, the scattering theory descrip- 
tion is noncausal with each interaction event equally as 
important as all others [6l|, [g^l • Note that this expression 
constitutes the Born approximation to the field pertur- 
bation as itself does not include any effects of field 
interaction. Dyson's expansion of subjected to the 
field would provide a systematic way to build in these 
effects perturbatively. 
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From Eq. ([70]) . the Born approximation to the inter- fields Eq and Ep is given by 
acting scattering T-matrix with respect to the external 



J2 W'j\{e^x ^ ■ (ple^r) 



rs—ij,ah\XY 

+ e(,+^ ■ (p|ex|r) 



s\eyi'\q) ■ Ey 
s|ex'|g) -E^^^llz/j) 



(+) 



(72) 



^ (^S|{e^"^ • {p\e^\r)[g%s{sl + fi^k)(.|ex'|g) • E 



(+) 



rs—ib,aj\XY 



e(.+) . (p|ex|r)[g«],,(eO - huj^,) {s\e^' \q) ■ 



where the off-diagonal matrix elements of the interact- 
ing molecular Green's function Qp in the last two terms 
were given in Eq. l|42|) . As discussed previously be- 
low Eq. the intermediate states associated with 
Raman scattering are restricted to particle-particle and 
hole-hole states; no particle-hole or hole-particle inter- 
mediate states contribute to its lowest-order theoreti- 
cal description. Consequently, in the following, we omit 
the off-diagonal components of Qp underlying the inter- 
acting scattering T-matrix. The first (diagonal) term 



stemming from hYip, (e^ -I- fiwk) represents the uncrossed 
scattering contributions analogous to those shown in 
Fig. [21 while the second (diagonal) term stemming from 
?il]*(e^ — Swk') represents the crossed terms (not shown) 
where a Raman-scattered photon is emitted in the ini- 
tial molecular state and an incident photon strikes in the 
final state. 

By including only particle-particle and hole-hole inter- 
mediate molecular-electronic states, the interacting scat- 
tering T-matrix reduces to 



/ / |/T-int(2) I \ 



E 

rs—ij.ab-^XY 

E 

rs—ij.,ab-^XY 



{p\ex^\r)5rs{s\ex"'\q) 



fil^k + £° - e° - mB]rr{^q + ^k) 

{p\ex^\r)5rs [n^)]rg(?it^k)(s|ea::"^|g) 



{p\ex''\r)5rs{s\ex'^\q) 



{p\ex'^\r)5rs\^f^)]rq{-fvjj\^i){s\ex'^\q) 



E ^ 

rs—ij^ab\XY 



^^^^''■''{i-[nfo)].,(?^k)?i[s*],^.(£g 



^k)}{i - [nfM-nw^')h[Y.l]PAel - r^c^k')} 



(+) 

(73) 



where two terms in the numerator which stem from find- 
ing a common denominator have been omitted in the 
third line. As with the off-diagonal terms, these terms 
involve higher powers of the molecular polarizability and 
correspond to repeated photon scattering events with the 
molecule. Identification of the inverse retarded HE po- 
larization propagator from Eq. (jB2[) in the denominator 
of this expression has been made and can be used to sim- 



plify the denominator in the last line as 

[nfo)],,(±;i^k,k')^P*]rr(eg ± ^k,k') 

= 5]Tr{5,^,„,(±ftc.k,k')-A 

tu 

de' ~ 
■ / — ft'5RPA(e')[Gfo)]tn(e° ± Tiw^m' + ' A}. 

(74) 



Within the RPA, the integral in this expression is the re- 
tarded component of the integral computed previously in 
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Eq. (ISS)) . Two different transition polarizabilities appear 
in Eqs. ^ and dTJ); in light of Eq. dBH) they are 



Performing the summation over X, Y results in the fol- 
lowing expression for the interacting scattering T-matrix 



pq^rs 



-n^k') = (p|ex|r)[n(^)]^<,(fiwk)(s|ex'|g) 
+ (p|ex|r)[n(^)],,(-;iwk')(s|ex'|g) 
5*^,„,(±fi^k,k') = (w|ex|r)[nfo)],,(±fiwk,k')(Hex'|t). 



; / |^int(2) I \ 



E 



rs=^,,ab {l - (fi^k)?i[S*]f, (eg + ;ic.k)}{l - [fifo) (-^c^k' ] (eg - /Jc^k')} 



e[, ^ . 5*^^,,(fii^k, -?iwk') • E^+^(cc;k) + E|,p^(-Wk') • fipl..(fi^k, -?iwk') • E^+^ 



r 

_ ^f'^J"'^ + 5^pA(-t^k')A'^^'5'"" + (^^^A"^ -^ 



{i-[n«)],,(Swk)^P 



*1P (pO 

RlrrK^q 



'OP (^k) 
'SRPA(^k) 



fii.k)}{l - [nf ],,(-fi^k')^P^]^.(e° - fi^k')} 



r 



(76) 



The numerator is of the form [1 + 5rpaA]^, while, in 
light of Eq. (jTi]) . the denominator is of the form [1 — 
Sa/ASrpaA]^; both are perfect squares when cjk' = ^k- 
Further, with Eq. (|74p . the following incident and Raman 
(anti-)Stokes shifted (unitless) enhancement factors [7, 
Il8l | may be rigorously defined by 



g'rq 



(^k) 



'^^^ + 4pA(-^k')A^-^ 



1 - [n«)],,(-fi^k')^[s^].^.(eg - 



1 - [nf^)U{rKo^)m*\Ue^^ + fia;k) 



(77) 



in terms of which the transition amplitude becomes 



X {v'.j\[(^MVplrr{^\^^ -ftWk')|i^,/)gr'g (^k)£^a 



(+) 



Due to this factorization, the quantum-mechanical 
Raman-scattering intensity associated with the interact- 
ing molecule-particle system displays a fourth-power en- 
hancement when both incident and Raman scattered 
fields share the same frequency, in analogy to the classi- 
cal case of two coupled dipoles discussed in Sec. |lTl see, 
e.g., Eq. Otherwise, g and g' do not maximally 

constructively multiply but, rather, contribute a factor 
of \gg'\'^ to the enhanced Raman-scattering intensity. 



Upon expanding the molecular-electronic polarizabil- 
ity OLM in Eq. (|78p around the molecule's equilibrium 
geometry Qo and keeping only the linear term, subse- 
quent application of Fermi's golden rule yields the follow- 
ing quantum-mechanical result for the SERS intensity 



:?SER^ ^ i^(A/-' + l)k-v| E4"^*(k')(^,7l(Qj - Qo) ■ VQ,[5M]2,V(^^k, -S^k')l^j)4^^"(k) 

lQ{K) C I j 



5iiSv<y + 4^PA(-c.k')A*^5''- + J^^A'^^ag.^AK) + 4^A(-^k')A^^A''^5g,;AK) 



[1 - EtnTr{<,„t(±u;k,k') • A • /(dc.'/2^z)5RPA(u;')[G« ]t„(eO/S ± a;k,k' + a;') • A}] 



(79) 



in the k' direction with polarization el^, ■* (k') from a single density fluctuations in a nearby classical metallic particle 
molecule that is interacting with its self-induced plasma 
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in the presence of an external radiation source propagat- 
ing in the k' direction with polarization el^^^(k). To our 
knowledge, this is the first place in the literature where a 
SERS intensity has been derived entirely from first princi- 
ples that explicitly treats the coupling and back-reaction 
effects of a quantum-mechanical molecular-electronic sys- 
tem with a nearby metallic particle supporting collective 
excitation of conduction electrons in the presence of a 
perturbing radiation field. Both image and local-field ef- 
fects, which contain the essential physics underlying the 
electromagnetic mechanism of SERS, have been general- 
ized beyond the classical model of Sec. |lT]to a quantum- 
mechanical setting. 

We note that it is not possible to exactly factor this 
general expression into a normal Raman-scattering com- 
ponent and an enhancement factor as was done in the 
classical case beyond what was already written in Eqs. 
([77)1 and ((78)) . This is due to the summation over the 
molecular states r and contraction on spatial (Greek) in- 
dices that connect all parts of the expression together. 
Nonetheless, it is still possible to identify a normal 
Raman-like component and an enhancement factor as the 
prefactor and quotient respectively in Eq. ([75)) . As in the 
case of normal Raman scattering, — Qo)|'^j) is 

nonzero whenever v'j — z/j ± 1. Note that, due to for- 
matting constraints, the square in the denominator is to 
be understood as the product of two separate factors: 
one associated with an incident photon with frequency 
and another associated with a Raman (anti-)Stokes 
scattered photon with frequency Wk' • 



E. Molecular and Particle Response 

In the previous sections we have developed both non- 
perturbative and perturbative expressions for the one- 
body molecular-electronic Green's function that built in 
the effects of interaction with a nearby metal particle un- 
der the influence of an external radiation field. Together 
with the scattering T-matrix it was demonstrated that 
these Green's function techniques underlie a quantum- 
mechanical picture of SERS that generalizes (and, in the 
appropriate limits, reduces to) the classical theory pre- 
sented in Sec. [TTl 

Until this point only perturbations acting upon the 
molecular system have been treated. Where the par- 
ticle's response to the molecule and field is prescribed 
by the RPA polarizability (|57p or some variant, e.g., a 
Drude polarizability, there is no need to explicitly quan- 
tify their perturbing effects. However, it is our desire to 
go beyond the RPA and solve for the classical dynamics 
of the particle's first-order induced dipole moment p*^^^ 
and resultant polarizability ap as set up by a nearby 
Raman-active molecule (both image and local-field ef- 
fects) and external radiation field. Incorporation of the 
latter interaction is straightforward, however, before in- 
cluding the former, we must first compute the induced 
molecular-electronic dipole moment d'-'^' of the interact- 



ing molecular system. Knowledge of d^^^ closes our theo- 
retical formulation of the many-body SERS problem and 
renders it computationally well-defined. It is our goal 
here to compute this molecular response and develop ba- 
sic governing equations for the dynamics of the metallic 
particle. 



1. Review of linear-response theory 

The linear response of a quantum many-body system 
[39I lisl . m, [ill subjected to a sufficiently weak external 
time-dependent perturbation V{t) ~ J p{x,t)v{:x.,t)cPx 
may be characterized by the first-order fiuctuations in 
its charge density p according to 

(/5(x,i))ext - P(x,t) =p(l)(x,^) + ••• 
= lId'^'H^dt'{[pi^\ t'), p(x, <)] )^;(x', t') -I- ■ • ■ , 

(80) 

where the expectation value (■ ■ • ) is taken within some 
unperturbed many-particle reference state, while the 
many-particle states underlying (• • • )cxt include the ef- 
fects of the external perturbation V; (• • • )oxt reduces to 
(• • • ) in the hmit of V ^ 0. 

For definiteness and continuity with the previous we 
henceforth choose the states underlying (• • • ) to be the 
noninteracting HP reference state |$Hp) and work in the 
interaction picture with respect to it. From Eq. (jSU)) . we 
see that the zeroth-order (static) unperturbed HP den- 
sity p — ('I'hfIpI'^'hf)' while the first-order density fluc- 
tuations 

(x, t) = i/d^x'nfo) (x, t; x', t')v{^', t') (81) 

can be reexpressed in terms of the retarded HP polariza- 
tion propagator 11^^ defined in Eq. (|B2j. Here, the time 
integral is extended to all times through the definition 

zHfo)(x, t; x', t') ^ 0{t - t')(1'HFl [P(x, t), p(x', t')] \Kf) 

= e{t-t'){^^p\[Sp{^,t),Spi^',t')]\KF), 

(82) 

where, as before, Sp = p — {^hf\p\^hf) = P ^ P- 

Taking the external perturbing potential f (x, t) — 
— (—ex) • Eo(x = 0,t) in the electric dipole-interaction 
approximation, we can derive an explicit expression for 
the first-order induced dipole moment d^^^ generated by 
V. It is 

d(i)(t) =/d3a;(-ex)p(i)(x,t) 

= ^/d3a;d4a;'(-ex)Hfo)(x,t;x',t')(-ex') •Eo(t') 
= Jdt'aM{t,t') -Eoit'), 

(83) 



21 



where the molecular-electronic HF polarizability (jB3|) 
may also be defined by 

~iha.M{t,t') = /d^a;d^x'(-ex)in(^)(x,t;x',t')(-ex') 
^e{t-t'){^U[5d{t),5d{t')\\^^^). 

(84) 



From Eqs. ([55]) and (|B2p . we see that, through the under- 
lying polarization propagator n^^, the polarizability cxm 
acts as a integral response kernel that depends upon the 
set of all single-particle electronic states of the molecule 
and parametrically upon all underlying nuclear coordi- 
nates. The external field Eq induces electronic density 
fiuctuations that oscillate with the field until the field 
frequency is close to an excitation energy of molecule. In 
this way, information on molecular electronically-excited 
states and transition amplitudes can be obtained, among 
other quantities. 



2. Molecular response to perturbations induced by a 
metallic particle and external electric field 

Equation ([5T|) prescribes a method to compute the 
first-order density fluctuations of a noninteracting molec- 
ular system induced by the external potential t'(x, t) = 
— (—ex) • Eo(t) from integration of v against the ker- 
nel n^^. The quantity p*^^-* describes how the molecular- 
electronic density alone responds to the perturbation v. 
However, we are not interested in knowing how the non- 
interacting molecular system responds to v, but, rather, 
want to study the response of an interacting molecule- 
particle system. Replacing the noninteracting integral 
kernel Il^-j with the interacting polarization propagator 

Hfj describing the coupling between molecular electrons 
and conduction electrons in the particle achieves pre- 
cisely this goal. In particular, it is desired for Ufj to 
incorporate those effects already built into the interact- 
ing molecular-electronic Green's function . 

A simple and general relationship exists between the 
polarization propagator and the one-body Green's func- 
tion. Application of Wick's theorem to the full Hm, de- 
fined in analogy to Eq. (jBl|) . reveals that 



«nM(x,t;x',i') 
= ($l^Fm^^(x, t)*(x, 0*^(x', t')*(x', i')}l*HF) 

-($ilF|*t(x,t)^(x,i)|$l^F)(*HFl*^(x',t')*(x',i')KF) 

(85) 

where the electron field operators are expanded onto the 
underlying basis of interacting orbitals Xg s-nd interact- 
ing electronic creation and annihilation operators Cq de- 
fined in Sec. IIIICll This result together with Eq. ((5T|l . 
demonstrates that the first-order density fluctuations of 



the interacting molecule-particle system may be com- 
puted from 



piy(x,0 - ■i/dVnf,(x,i;x',Ot'(x',0 

= ^/^'^' [^''(^' ^Oe'^Cx', t'; X, t)] ^v{x', t'), 

(86) 

or, alternatively, from the spectral form 
= i/dVn« (x,x';^)?J(x',u;) 
= ^Jd^x' [/^^''(x, x'; CO + iu')gP{^', x; cj')] ^u(x', c.) 



1 



(87) 



with pg-matrix elements 



^E[/^^"^(-+-')ei;(-')]/.«M 



From this last expression it is clear that knowledge of 
may be attained once the (retarded component of 
the) interacting molecular polarization propagator 11^ 
is known. By assuming that the self energy varies 
so slowly with frequency that it is appropriate to ap- 
proximate [E*]pp(u;) ^ [T,p]pp{e°/h) and [S*]pg(w) w 
[Y.p]pq{ep/h), the interacting 11^^ reduces to 



dpfSgq 



Lo-el/h + e^Jh- 



A\qq 



SprSgq 



- el/h + el/h - [A^p - (^/2)rg -I- [A^, -f (^/2)r^J 

(89) 

in the weak-coupling limit. The shifting (A^) and broad- 
ening (r-'^) of the single-particle states p and s due to 
the explicit treatment of interaction have been quanti- 
fied previously in Eqs. ((60)l and (|6T|) . 

Together with Eq. ([55]) . this approximate expres- 
sion for the interacting molecular-electronic polarization 
propagator can be used to compute the first-order elec- 
tronic density fluctuations in the molecule induced 
by the external fleld Eq and, through , by interac- 
tion with molecule-induced density fluctuations um in a 
nearby classical metal particle. In the dipole approxima- 
tion, the electric fleld at the point X2 generated from p^^-* 
located at the origin is given by 



Em(x2,c^) = A-di^ H = A./(-ex)p\V(x,c^)d'x (90) 

in the frequency domain. It will be now be demonstrated 
how the induced local molecular electric fleld Em in- 
fluences the dynamics of the conduction electrons in a 
nearby particle. 
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3. Particle response to perturbations induced by molecular 
electrons and external electric field 

In order to expose certain notable features of our for- 
malism, such as the enhanced Raman-scattering intensity 
in Eq. (j79p . we have imposed a predetermined dynamics 
upon the collective excitations of conduction electrons in 
the metal particle: namely, that of the RPA high-density 
electron gas ([571) . This ad hoc choice has facilitated ana- 
lytical derivation. Here we provide the minimal theoret- 
ical framework necessary to go beyond such a formalism 
and treat the dynamics of the metal particle's conduction 
electrons explicitly. 

In Eq. (fTO|) of Sec. IIIIl the many-body Hamilto- 
nian Hp of the particle and its interaction with a nearby 
Raman-active molecule under the influence of an exter- 
nal electric field was presented. Its external interaction 
component can be decomposed according to 

HLit) « -P(t) • [Eo(x2,t) + A • d«(t)], (91) 

where the interaction potential W ~ — p- A-d^^ has been 
multipole expanded and only the term of dipole order was 
retained. Linear-response theory may now be employed 
to compute the electronic density fluctuations induced in 
the particle by interactions with the external field Eq and 
with the dipole field of the molecule A ■ d^^ . This latter 
term represents the local- field effect of the molecule upon 
the particle. 

By discretizing the particle's associated electric dipole 
fluctuations onto a three-dimensional spatial grid, the re- 
sulting response equation for the jth induced dipole mo- 
ment is given by 

p^^\t)^ Jdt'ap,it,t') 

•[Eo,(0 + E^^fe-pi'^(^') + A,-diy(t')], 

(92) 

where Aj«f = [3r|fcfJ^ - (5«'"]/r|fc with rjk = |x2j - X2fc| 
is the dipole tensor associated with the interaction of 
two such discretized induced dipoles in the particle la- 
beled by j,k {k ^ j), and Af = [3f|fJ - d^^j/r^ with 
Tj = |xi — X2j | is the dipole tensor associated with the in- 
teraction of the discretized induced dipole in the particle 
labeled by j with the induced molecular dipole moment 
d^-* located at Xi, which may be taken as the coordi- 
nate origin. The underlying particle polarizability apj 
in Eq. ([5^ may be approximated by an appropriately 
discretized polarizability that is consistent with the op- 
tical theorem [6^. Both image and local-field effects are 

included in this expression by the terms J2k^j ^'jk ' Pfc^^ 
and Aj • d^^'' respectively. 



In Fourier space, p^- takes the form 

p«(a;) = 5p,(c^)-[Eo,(c^)+5] A;.,.p«(c^)-HA,.diV(a;)], 

(93) 

which may be inverted to yield 

[ap]{u;)S,k-{l-5,k)A'^^]pi'\io) = Eo,{u;)+EMj{io). 

k 

(94) 

Hence, discretization leads to a linear AX — B system 
of equations that are solvable by standard linear algebra 
routines [6^. Similar response equations are well known 
in the literature; see, e.g., Ref. [6^. This spatial dis- 
cretization enables the practical treatment of anisotropic 
metal particles of arbitrary shape and size. 

In the dipole approximation, the (continuum) electric 
field stemming from the solution of these equations may 
be written as Ep = A-p^^^^. By allowing the induced elec- 
tric field of the particle's conduction electrons to act back 
upon the molecular system through, e.g., Eq. (jSOp . the 
many-body formalism presented in this article is math- 
ematically closed and well defined. Iteration between a 
quantum-chemical molecular response calculation where 
the incident field Eq and induced field of the particle Ep 
enter as perturbations and a classical-electrodynamical 
metal particle response calculation where the incident 
field Eq and induced field of the molecule Eji/ enter as 
perturbations should be performed until self consistency 
is reached between the two systems. 

IV. SUMMARY AND CONCLUSIONS 

We have presented a unified and didactic approach to 
the understanding of the microscopic theory of single- 
molecule SERS from a nanoscale metal particle at zero 
temperature. Nonperturbative and perturbative many- 
body Green's function techniques are employed to build 
in the interaction between a molecular-electronic system 
and the conduction electrons of a nearby metallic par- 
ticle in the presence of an external radiation field from 
first principles. Both image and local-field effects be- 
tween molecule and metal are explicitly included. Due to 
its generality, other relevant approaches from the litera- 
ture, including those that are purely classical and purely 
quantum-mechanical, are obtained by taking appropriate 
limits of our formalism. 

With emphasis placed upon practical initial numer- 
ical implementation, molecular-electronic correlation is 
restricted to the level of Hartree-Fock mean-field theory, 
while a Bogoliubov decomposition of the metallic parti- 
cle's plasmon field is effected to reduce these collective 
electronic excitations to a classical field; extension of for- 
mer through a Kohn-Sham density-functional or M0ller- 
Plesset perturbation theory is discussed; specialization 
of the latter to the RPA allows for the analytic presenta- 
tion of several salient features of the theory such as the 



23 



enhanced Raman-scattering intensity in Eq. ((79)) and in- 
cident and Raman (anti-)Stokes scattered enhancement 
factors in Eq. ([77j) . We believe that such expressions, 
where interaction effects are systematically included from 
first-principles, have never before appeared in the litera- 
ture. 

Going beyond the RPA, we explicitly describe the re- 
sponse of the metal particle's conduction electrons to the 
external field and to the induced electronic density in 
a nearby molecule by discretizing the associated classi- 
cal response equations onto a spatial grid. By iteration 
until self consistency is reached, the reaction and back- 
reaction effects between molecular and particle systems 
with each other and with the external perturbing radia- 
tion field may be incorporated: thereby mathematically 
closing the theory. Having laid out our many-body for- 
malism in this article and demonstrated its relation to 
and generalization of other approaches from the litera- 
ture, as well as some of its important properties, imple- 
mentation is currently underway to numerically realize 
these equations and provide theoretical support to cur- 
rent single-molecule SERS experiments. 



APPENDIX A: INCLUSION OF 
ELECTRON-ELECTRON INTERACTION 
EFFECTS 

The effects of electron-electron interaction stemming 
from Hint in Eq. (jl9p may be straightforwardly included 
in the above formalism with either M0ller-Plesset per- 
turbation theory or density-functional theory. In the lat- 
ter case, the HE orbital equation (PT|) is replaced by the 



J 



zn(«)(x,i;x',0 = {^^p\T{dp{^,t)dpi^',t')}\Kp) 

= J2 0;(x)0p(x)inW„(t,t')0;(x')0.(x') 

pqrs 



Kohn-Sham orbital equation 

0- [;io(x) + C/o(x)-ef ](^,(x) 

+ ^(fc|y|fc)(x)^,(x) + (5i?^^M/5p(x))^,(x), 
k 

(AI) 

where the total molecular-electronic energy 



It has the retarded component 



Em = / p(x)[C/o(x) + VeA^)]d^x 

+ {l/2)Jp{x.)V{x,x')p{x')d^xd^x' + G[p] 



(A2) 



may be expressed as the sum of the external electron- 
nuclear attraction V^xt and universal functional 



G[p]^To[p]+E^c[p] 



(A3) 



with noninteracting electronic kinetic energy To = To [p] . 
If, in addition, E^fj included the interaction effects 
among conduction electrons within the metal particle (or 
a subset thereof) and between these conduction electrons 
and the molecular-electronic system, then the present ap- 
proach (with a few minor modifications) would recover 
the work presented in Refs. [H [H, [Ij, HI]. In this 
way, the chemical mechanism of SERS can additionally 
be incorporated within the present formalism. 
APPENDIX B: POLARIZATION PROPAGATOR 
AND LINEAR POLARIZABILITY 

The noninteracting HE polarization propagator (ssl . 
I39I I is defined by the expectation 



(1 - P>' 



u; + iO+ + (eo - eO)/h lu - iO+ + (e^ - eO)/h 



}0;(x')0.(x'). 
(Bl) 



mfo)(x, i; x', t') = 9{t - t')($HFl ['^/5(x, t), 6p{^, t')] |$^^p) 

= ^ 0;(x)^p(x)z[n«)]p,,,(t,t')'/';(x')0.(x') 

pqrs 

. duj 



(B2) 



pqrs 
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expressed in terms of the commutator [■,■]. It represents 
the probabihty amphtude that a molecular-electronic 
density disturbance 5p originating at the space-time 
point (x',t') will be found later within the iV-electron 
Fermi vacuum |<i>Hp) at the space-time point (x, t), where 
the electron density operator /5(x, i) = >I'^(x. i)>I'(x, <) is 
decomposed into static and interaction-induced excita- 
tion parts according to p = ('I'hfI/^I'^hf) + ^P- Pi'oof of 
Eq. (|3ip may be demonstrated by retaining the plasmon 
field operators £7 to second order in perturbation the- 
ory in Eq. (I30|) . defining plasmon Green's functions, and 
then effecting the Bogoliubov decomposition ^ together 
with the Heaviside identity 6{x) -t- 0{—x) = 1. 

From Eq. (jB2p . the linear HF polarizability is defined 

by 

-ihai^{t,t') = - exi\pH^io)Urs{t,t')(r\ - ex'^ls) 

pqrs 

(g| - ex"'\p)(j)\ - ex^\q) - 

(B3) 

with interaction-induced dipole moment 5d{t) — 
J {—ex.)Sp{x,t)(Px being related to its expectation value 
in an externally perturbed reference state (■ • • )oxt by 
(M(t))ext = d(i) (t) + • • • = Jdt'aMit,f) ■ Eo(t') + ■ • • ; 
see Sec. IIII E II for a review of linear-response theory, 
and, in particular Eq. (j83p . It can be shown that aM 
satisfies [Sm(w)]* = cj). We note that there has 

been an ongoing debate in the literature regarding the 
damping sign (±iO^) convention in the linear (Kramers- 
Heisenberg) polarizability. This issue has recently been 
positively resolved by Bialynicki-Birula and Sowihski [Blj 
for linear processes and by Mukamel ^62] for nonlinear 
processes, and the differences between scattering and re- 
sponse points of view have been clarified. To briefiy 
summarize, both signs are correct but apply to differ- 
ent physical situations [6l|. The fundamental principle 
of causality is always respected, however, response theory 
describes interaction processes occurring with a particu- 
lar causal time ordering where the initial state and form 
of interaction are specified and a sum over final states 



is performed, while scattering theory equally includes all 
possible time orderings of interactions between specific 
initial and final states and, in this sense, is noncausal 
(6^ . The response functions [e.g., Eq. (jB3|) ] and tran- 
sition amplitudes [e.g., Eq. (jB4p ] computed here from a 
many-body perspective are consistent with Refs. [6ll[63|: 
the differences in damping sign between Eq. (|B3|) and Eq. 
(IB4[) below are correct. 

In order to compute Raman transition amplitudes, we 
introduce the Lehmann representation of the generalized 
nonintcracting molecular transition polarizability 

= X] [^^1 ~ ex^\r)Srs[fifo)]rqii^){s[ - ex'" \q) 

rs 

+ {p[ - ex"'\r)5rs[^fQ)\rq{~uj'){s\ - ex^\q) 

_ r (p| - ex^\r)5rs{s[ - ex"^\q) 
-^l (£°-eO)A + c. + zO+ 

{p[-ex"'[r)drs{s[-ex^q) - 

(B4) 

between the single-particle states q and p with two po- 
tentially different frequencies u> and w', defined in terms 
of the rg- matrix elements of the retarded HF polarization 
propagator as 

[^fo)UpsH = SrpS^M - p'j.mfo^Uiu;). (B5) 

We point out that the sums in all previous expressions 
for the polarization propagator and polarizability extend 
over the set of all single-particle states of the molecule 
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